Smirnov's fermionic observable away from criticality 



V. Beffara H. Duminil-Copin 

June 29, 2011 

Abstract 

In a recent and celebrated article [TBI , Smirnov defines an observable for the self-dual random- 
cluster model with cluster weight q = 2 on the square lattice Z 2 , and uses it to obtain conformal 
invariance in the scaling limit. We study this observable away from the self-dual point. From 
this, we obtain a new derivation of the fact that the self-dual and critical points coincide, which 
implies that the critical inverse temperature of the Ising model equals \ log(l + \2). Moreover, 
we relate the correlation length of the model to the large deviation behavior of a certain massive 
random walk (thus confirming an observation by Messikh |14|). which allows us to compute it 
explicitly. 

Introduction 

The Ising model was introduced by Lenz |12j as a model for ferromagnetism. His student Ising 
proved in his PhD thesis [9] that the model does not exhibit any phase transition in one dimension. 
On the square lattice L = (Z 2 , E) , the Ising model is the first model where phase transition and 
non-mean- field behavior have been established (this was done by Peierls |15]). 

An Ising configuration is a random assignment of spins { — 1, 1} on Z 2 such that the probability 
of a configuration a is proportional to exp[— (3 J2 a ~b a ( a ) a (b)], where j3 is the inverse temperature 
of the model and a ~ b means that (a, b) is an edge of the lattice, i.e. (a,b) £ E. Kramers 
and Wannier identified (without proof) the critical temperature where a phase transition occurs, 
separating an ordered from a disordered phase, using planar duality. In 1944, Kaufman and Onsager 
|10j computed the free energy of the model, paving the way to an analytic derivation of its critical 
temperature. In 1987, Aizenman, Barsky and Fernandez [TJ found a computation of the critical 
temperature based on differential inequalities. Both strategies are quite involved, and the first goal 
of this paper is to propose a new method, relying only on Smirnov's observable: 

Theorem 1. The critical inverse temperature of the Ising model on the square lattice Z 2 is equal 
to 

p c = i In (l + v^) ■ 

Beyond the determination of the critical inverse temperature, physicists and mathematicians 
were interested in estimates for the correlation between two spins, E ( g[o"(a)a'(6)] (where E^ denotes 
the Ising measure). McCoy and Wu [13] derived a closed formula for the two-point function, and 
an asymptotic analysis shows that it decays exponentially fast when f3 < f3 c . In addition to this, 
it was noticed by Messikh [TJ] that the rate of decay is connected to large deviations estimates for 
the simple random walk. In this article, we present a direct derivation of this link, which provides 
a quick proof of the following theorem: 
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Theorem 2. Let (3 < (3 C and let denote the (unique) infinite-volume Ising measure at inverse 
temperature f3; fix a = (ai, 02) £ L. Then, 

lim In (Efl[a(0W(na)]) = 01 arcsinhsai + a2 arcsinhsa2, 

n— >oo 7i 

where s solves the equation 

sf\ + (sai) 2 + v 7 ! + (sa2) 2 = sinh 2/3 + sinh" 1 2,8. 

Instead of working with the Ising model, we rather deal with its random- cluster representation 
(known as the random- cluster model with cluster weight q = 2). It is well known that one can couple 
this model with the Ising model (see e.g. [8] for a comprehensive study of random-cluster models). 
The spin correlations of the Ising model get rephrased as cluster connectivity properties of their 
random-cluster representations, which allows for the use of geometric techniques. For instance, the 
determination of j3 c is equivalent to the determination of the critical point p c for the random-cluster 
model. 

The understanding of the two-dimensional random-cluster model with q = 2 has recently pro- 
gressed greatly Qj]13]- thanks to the use of the so-called fermionic observable introduced by Smirnov 
|16j . which was instrumental in the proof of conformal invariance. This observable is defined on the 
edges of a finite domain with Dobrushin boundary conditions (mixed free and wired, see Section 1 
for a formal definition), and it is discrete holomorphic at the self-dual point p sc [ = \/2/(l + y/2). 

The idea of our argument is the following. Below the self-dual point, the observable can still 
be defined but discrete holomorphicity fails and the observable decays exponentially fast in the 
distance to the wired boundary. Along the free boundary, the modulus of the observable can be 
written exactly as a connection probability, so in the p < p s d regime the two-point function is 
exponentially small as well, and that implies that the system is then in the subcritical regime, thus 
providing the lower bound p c ^ p s d on the critical parameter. Theorem [T] then follows from duality. 
In fact, the rate of exponential decay (and therefore Theorem[2]) can be derived by comparing the 



observable to the Green function of a massive random walk (Proposition 4.1 ); the key ingredient is 
the observation that the observable is massive harmonic in the bulk for p < p S( i- The correspondence 
between the two-point function of the Ising model and that of the massive random walk was 
previously noticed by Messikh |14j . 

In Section [TJ we remind the reader of a few classic features of the random-cluster model. In 
Section[2j we define this observable away from criticality and gather some of its important properties 
— for instance, the fact that the observable on a graph is related to connection properties for sites 
on the boundary. In Section [3j we derive Theorem [T] by showing that the observable decays 
exponentially fast. Section [4] is devoted to a refinement of estimates on the observable, which leads 
to the proof of Theorem [2} 



1 Basic features of the model 



The Ising model on the square lattice admits a classical representation through the so-called 
random- cluster model with q = 2. This model can be studied using geometric arguments which are 
classic in the theory of lattice models. We list here a few basic features of random-cluster mod- 
els; a more exhaustive treatment (together with the proofs of all our statements) can be found in 
Grimmett's monograph [8j. Readers familiar with the subject can skip directly to the next section. 
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Definition of the random-cluster model. The random-cluster measure can be defined on any 
graph. However, we will restrict ourselves to the square lattice, denoted by L = (Z 2 ,E) with Z 2 
denoting the set of sites and E the set of bonds. In this paper, G will always denote a connected 
subgraph of L, i.e. a subset of vertices of Z 2 together with all the bonds between them. We denote 
by dG the (inner) boundary of G, i.e. the set of sites of G linked by a bond to a site of Z 2 \ G. 

A configuration uj on G is a random subgraph of G, having the same sites and a subset of its 
bonds. We will call the bonds belonging to to open, the others closed. Two sites a and b are said to 
be connected (denoted by a -<->• b), if there is an open path — a path composed of open bonds only 
— connecting them. The (maximal) connected components will be called clusters. More generally, 
we extend this definition and notation to sets in a straightforward way. 

A boundary condition £ is a partition of dG. We denote by to U £ the graph obtained from the 
configuration uj by identifying (or wiring) the vertices in £ that belong to the same class of £. A 
boundary condition encodes the way in which sites are connected outside of G. Alternatively, one 
can see it as a collection of abstract bonds connecting the vertices in each of the classes to each 
other. We still denote by wU( the graph obtained by adding the new bonds in £ to the configuration 
to, since this will not lead to confusion. Let o(to) (resp. c(w)) denote the number of open (resp. 
closed) bonds of uj and k(uo, £) the number of connected components of uj U £. The probability 
measure 4>\ q c °f the random-cluster model on a finite subgraph G with parameters p E [0, 1] and 
q £ (0, oo) and boundary condition £ is defined by 
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for any subgraph uj of G, where Z^ G is a normalizing constant known as the partition function. 
When there is no possible confusion, we will drop the reference to parameters in the notation. 

The domain Markov property. One can encode, using an appropriate boundary condition £, 
the influence of the configuration outside a sub- graph on the measure within it. Consider a graph 
G = (V, E) and a random-cluster measure (fit G on it. For F C E, consider G' with F as the set 
of edges and the endpoints of it as the set of sites. Then, the restriction to G' of (ffi G conditioned 
to match some configuration uj outside G' is exactly 4>p qG n where £ describes the connections 
inherited from uj U tp (two sites are wired if they are connected by a path in uj U tp outside G' - 
see (4.13) in [8]). This property is the direct analog of the DLR conditions for spin systems. 

Comparison of boundary conditions when q ^ 1. An event is called increasing if it is 
preserved by addition of open edges. When q ^ 1, the model is positively correlated (see (4.14) in 
[8]), which has the following consequence: for any boundary conditions ip ^ £ (meaning that ip is 
finer than £, or in other words, that there are fewer connections in ip than in £), we have 

< q d A ) < <4,g04) (1-2) 

for any increasing event A. This last property, combined with the Domain Markov property, 
provides a powerful tool in order to study how events decorrelate. 

Examples of boundary conditions: free, wired and Dobrushin. Two boundary conditions 
play a special role in the study of random-cluster models: the wired boundary condition, denoted 
by (pi, G , is specified by the fact that all the vertices on the boundary are pairwise connected; the 
free boundary condition, denoted by G , is specified by the absence of wirings between boundary 
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sites. These boundary conditions are extremal for stochastic ordering, since any boundary condition 
is smaller (resp. greater) than the wired (resp. free) boundary condition. 

Another example of boundary condition will be very useful in this paper. The following def- 
inition is deliberately not as general as would be possible, in order to limit the introduction of 
notation. Let G be a finite subgraph of L; assume that its boundary is a self-avoiding polygon in L, 
and let a and b be two sites of dG. The triple (G, a, b) is called a Dobrushin domain. Orienting its 
boundary counterclockwise defines two oriented boundary arcs ab and ba; the Dobrushin boundary 
condition is defined to be free on ab (there are no wirings between boundary sites) and wired on 
ba (all the boundary sites are pairwise connected). We will refer to those arcs as the free arc and 
the wired arc, respectively. The measure associated to this boundary condition will be denoted by 

<Pp,q,G or sm W 9g ■ 

Planar duality for Dobrushin domains. One can associate to any random-cluster measure 
with parameters p and gona Dobrushin domain (G, a, b) a dual measure. First, define the dual 
graph G* as follows: place a site in the center of every face of G and every face of L adjacent to 
the free arc, see Figure [TJ Bonds of the dual graph correspond to bonds of the primal graph and 
link nearest neighbors. Construct a bond model on G* by declaring any bond of the dual graph to 
be open (resp. closed) if the corresponding bond of the primal lattice is closed (resp. open) for the 
initial random-cluster model. The new model on the dual graph is then a random-cluster measure 
with parameters p* = p* (p, q) and q* = q satisfying 

(1 — p)q p*p 

P*{p,q) ■= 7^ v — : — , or equivalents — r = q, 

(l-p)q+p (l-p*)(l-p) 

with wired boundary condition on the dual arc adjacent to ab, and free boundary condition on the 
dual arc adjacent to ba. In particular, it is again a random-cluster model with Dobrushin boundary 
condition. This relation is known as planar duality. It is then natural to define the self-dual point 
Psd = Psd(q) by solving the equation p*(p sd , q) = p sd , which gives 

Psd(q) ■= z ^p . 

i + Vq 

This notion of duality has a natural counterpart, with the same formal definition, for free 
boundary conditions: the dual model is then a random-cluster model with parameters p* and q, 
with wired boundary condition. 



Infinite- volume measures and the critical point. The domain Markov property and com- 
parison between boundary conditions allow to define infinite-volume measures. Indeed, one can 
consider a sequence of measures on boxes of increasing sizes with free boundary conditions. This 
sequence is increasing in the sense of stochastic domination, which implies that it converges weakly 
to a limiting measure, called the random-cluster measure on L with free boundary conditions (de- 
noted by <j)p g)- This classic construction can be performed with many other sequences of measures, 
defining a priori different infinite-volume measures on L. For instance, one can define the random- 
cluster measure with wired boundary conditions, by considering the decreasing sequence of 
random-cluster measures on finite boxes with wired boundary condition. 

For our purpose, the following example of infinite- volume measure will be important: we define a 
measure on the strip Si = Zx [0, £]. The sequence of measures (0r_^ ^m}x[o°i]^ m ^ 1S increasing, in the 
sense that for any cylindrical increasing event A defined in the strip, the sequence (<f>\™m]^o£i(^)) 
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is well-defined for m large enough and is non-decreasing. This implies that the sequence of measures 
converges weakly as m goes to infinity. The limit is called the random-cluster measure on the infinite 
strip with free boundary conditions on the top and wired boundary condition on the bottom, and 
we will denote it by <jy£' 00 . 

When defining such measures in infinite volume by thermodynamical limits, it is natural to ask 
whether the limit depends on the choice of domains and boundary conditions used to build it; in 
the case of the random-cluster model, a more specific version of the question is whether taking free 
or wired boundary conditions affects the limit — these two being extremal, if the limits match, 
this implies uniqueness of the infinite-volume limit for all boundary conditions. It can be shown 
that for fixed q ^ 1, uniqueness can fail only on a countable set T> q of values of p, see Theorem 
(4.60) of [8j. From that (or rather from the weaker statement that the set of values of p at which 
uniqueness holds is everywhere dense in [0,1]), and from the fact that measures for larger values 
of p dominate those for smaller values, it is not difficult to show that there exists a critical point 
p c such that for any infinite- volume measure with p < p c (resp. p > p c ), there is almost surely no 
infinite component of connected sites (resp. at least one infinite component). Moreover, it is also 
known that the infinite- volume measure is unique when p < p s ^. 

Remark 1.1. Physically, it is natural to conjecture that the critical point satisfies p c = p s d- Indeed, 
if one assumes p c 7^ p s d, there should be a phase transition due to the change of behavior in the 
primal model at p c and a second (different) phase transition due to the change of behavior in the 
dual model at p*. This is unlikely to happen — in fact, constructing a natural- looking model 
exhibiting two phase transitions is not so easy; but the equality of p c and p s d is only known to hold 
in a few specific cases. 

In the case of the random-cluster model on the square lattice, the authors proved recently [2] 
that indeed p c (q) = Psdil) f° r all g ^ 1 (therefore determining the critical temperature for all 
(/-state Potts models on L). The argument does not use Smirnov's observable, but it is quite a bit 
longer than the one we present here, is not as self-contained (mostly because it depends on recent 
sharp-threshold results by Graham and Grimmett [HI [7] ) , and it provides less information on the 
subcritical phase. 

Coupling with the Ising model. The random-cluster model on G with parameter q = 2 is of 
particular interest since it can be coupled with the Ising model; consider a configuration to sampled 
with probability (f)® 2 G and assign independently a spin +1 or —1 to every cluster with probability 
1/2. We are now facing a model of spins on sites of G. It can be proved that the law of the 
configuration corresponds to the Ising model at temperature (5 = (3{p) = — |ln(l — p) with free 
boundary condition. 

We are then equipped with a "dictionary" between the properties of the random-cluster model 
with q = 2 and those of the Ising model. One instance of this relation is given by the useful identity 



where the left-hand term denotes the correlation between sites and a for the Ising model at inverse 
temperature (3 on the graph G with free boundary condition. 

The critical inverse temperature j3 c of the Ising model is characterized by the fact that the two- 
point correlation undergoes a phase transition in its asymptotic behavior: below /3 C , the correlation 
goes to when a goes to infinity, while above it, it stays bounded away from 0. The previous 
definition readily implies that (3 C = — 2log(l — p c (2)). In order to prove Theorem [IJ it is thus 
sufficient to determine p c (2). Notice that the inverse temperature corresponding to the self-dual 




(1.3) 



5 



point is given by (3(p s d) = \ ln(l + y/2) so that what needs to be proved can be written as p c (2) = 
Psd{2). 

The same reasoning implies that we can compute correlation lengths for the random-cluster 
model in order to prove Theorem [2| 

2 Definition of the observable 

From now on, we consider only random-cluster models on the two-dimensional square 
lattice with parameter q = 2 (we drop the dependency on q in the notation). 

The medial lattice and the loop representation. Let (G, a, b) be a Dobrushin domain. In 
this paragraph, we aim for the construction of the loop representation of the random-cluster model, 
defined on the so-called medial graph. In order to do that, consider G together with its dual G*; 
declare black the sites of G and white the sites of G*. Replace every site with a colored diamond, 
as in Figure [!} The medial graph G<> = (V<>, E ) is defined as follows (see Figure [T] again) : E is 
the set of diamond sides which belong to both a black and a white diamond; V<> is the set of all 
the endpoints of the edges in E . We obtain a subgraph of a rotated (and rescaled) version of the 
usual square lattice. We give G an additional structure as an oriented graph by orienting its edges 
clockwise around white faces. 
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Figure 1: Left: A graph G with its dual G* . The black (resp. white) sites are the sites of 
G (resp. G*). The open bonds of G (resp. G*) are represented by solid (resp. dashed) black 
bonds. Right: Construction of the medial lattice and the loop representation: the loops are 
interfaces between primal and dual clusters. 

The random-cluster measure with Dobrushin boundary condition has a rather convenient rep- 
resentation in this setting. Consider a configuration u), it defines clusters in G and dual clusters 
in G* . Through every vertex of the medial graph passes either an open bond of G? or a dual open 
bond of 67*. Hence, there is a unique way to draw Eulerian (i.e. using every edge of E exactly 
once) loops on the medial lattice such that the loops are the interfaces separating primal clusters 
from dual clusters. Namely, a loop arriving at a vertex of the medial lattice always makes a ±7r/2 
turn so as not to cross the open or dual open bond through this vertex (see Figure [l] yet again). 

Besides loops, the configuration will have a single curve joining the vertices adjacent to a and 
b, which are the only vertices in V with three adjacent edges within the domain (the fourth edge 
emanating from a, resp. b will be denoted by e a , resp. eb). This curve is called the exploration path; 
we will denote it by 7. It corresponds to the interface between the cluster connected to the wired 
arc and the dual cluster connected to the free arc. 



6 



This gives a bijection between random-cluster configurations on G and Eulerian loop configu- 
rations on G<>. The probability measure can be nicely rewritten (using Euler's formula) in terms of 
the loop picture: 



\oj) = = , where x(p) : 



Z(p,G) ' (1~P)V2 



and Z(p, G) is a normalizing constant. Notice that p = p s d if and only if x(p) = 1. This bijection 
is called the loop representation of the random-cluster model. The orientation of the medial graph 
gives a natural orientation to the interfaces in the loop representation. 

The edge observable for Dobrushin domains. Fix a Dobrushin domain (G,a,b). Follow- 
ing [16] . we now define an observable F on the edges of its medial graph, i.e. a function F : E — > C. 
Roughly speaking, F is a modification of the probability that the exploration path passes through 
an edge. 

First, introduce the following definition: the winding Wr(z, z') of a curve V between two edges 
z and z' of the medial graph is the total rotation (in radians and oriented counter-clockwise) that 
the curve makes from the mid-point of edge z to that of edge z', see Figure [2} We define the 
observable F for any edge e G E as 

F(e):=^ 6 fe^( e ^l ee7 ), (2.1) 



G y '^eG 7 , 

where 7 is the exploration path. 

Remark 2.1. In [16J, Smirnov extends the observable to vertices — as being the sum of F on 
adjacent edges — in order to study the critical regime. Properly rescaled, this function converges 
to a holomorphic function, which is a key step towards the proof of conformal invariance. Away 
from criticality, it is more convenient to work directly with the observable on edges. 

The following three lemmas present the properties of the observable we will be using in the 
proofs of both theorems. They have direct counterparts in Smirnov's article [16] (in particular, the 



idea of the proof of Lemma 2.2 can be found in the proof of Lemma 4.12 of [IE]), and as such they 
are not completely new. We still include their proofs here since our goal is to keep the present 
paper as self-contained as possible. 

Lemma 2.2. Let u G G be a site on the free arc, and e be a side of the black diamond associated 
to u which borders a white diamond of the free arc, see Figure^ Then, 

\F(e)\ = 4> a (j b {u o wired arc). (2.2) 



Proof. Let u be a site of the free arc and recall that the exploration path is the interface between 
the open cluster connected to the wired arc and the dual open cluster connected to the free arc. 
Since u belongs to the free arc, u is connected to the wired arc if and only if e is on the exploration 
path, so that 

(f>Q b (u wired arc) = 4>Q b {e G 7). 

The edge e being on the boundary, the exploration path cannot wind around it, so that the winding 
(denoted W\) of the curve is deterministic (and easy to write in terms of that of the boundary itself). 
We deduce from this remark that 

\F{e)\ = |^(e^ w ^l ee7 )| = M w ^ 6 (e G 7 )| 

= <j>Q b {e G 7) = <f>Q b (u O wired arc). □ 
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Figure 2: Left: A schematic picture of the exploration path and a boundary point u, together 
with two possible choices ei and e-z for e. If u is connected to the wired arc, the exploration 
path must go through e. Right: The winding of a curve. In the first example, the curve did 
one quarter-turn on the left and one quarter-turn on the right. 



For a random-cluster model, one can use the parameters p or x interchangeably. We introduce 
a third parameter which will be convenient: let a = a(p) £ [0, 2tt) be given by the relation 

^ ■= £rfri- <M) 

e™i^x{p) + 1 

Observe that a{p) = if and only if p = p s d and a{p) > for p < p s d. With this definition: 

Lemma 2.3. Consider a vertex v £ K> with four adjacent edges in E . For every p 6 [0, 1], 

F(A) + F{C) = e ia(p) [F{B) + F{D)\ (2.4) 

where A and C (resp. B and D) are the adjacent edges pointing towards (resp. away from) v, as 
depicted in Figure^ 

Proof. Let v be a vertex of K> with four adjacent edges, indexed as mention above. Edges A and 
C play symmetric roles, so that we can further require the indexation to be in clockwise order (see 
one such indexation in Figure [3]) . Recall that any vertex in V corresponds to a bond of the primal 
graph and a bond of the dual graph. We consider the involution s on the space of configurations 
which switches the state (open or closed) of the bond of the primal lattice corresponding to v. 

Let e be an edge of the medial graph and denote by e w := (j)Q b (uj) e2 W T'( e > e f>) l eg7 the contribution 
of uj to -F(e). Since s is an involution, the following relation holds: 



:= ^ = \ \- eu) + e 



«(w)J 



In order to prove (2.4), it suffices to prove the following for any configuration uj: 

Aj + A a(u) + C u + C 9{u) = e ia ( p ) [B u + B a(u) +D U + D s(u}) ) . (2.5) 
When 7(0; ) does not go through any of the edges adjacent to v, it is easy to see that neither does 



j(s(uj)). All the contributions then vanish and (2.5) trivially holds. Thus we can assume that j(uj) 



passes through at least one edge adjacent to v. The interface follows the orientation of the medial 
graph, and thus can enter v through either A or C and leave through B or D. Without loss of 
generality we assume that it enters first through the edge A and leaves last through the edge D; 
the other cases are treated similarly. 
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Figure 3: Left: Indexation of the edges adjacent to v. Right: Two associated configurations 
uj and s(w). In this picture, v corresponds to a vertical bond of the primal lattice. 



Two cases can occur: either the exploration curve, after arriving through A, leaves through B 
and then returns a second time through C, leaving through D; or the exploration curve arrives 
through A and leaves through D, with B and C belonging to a loop. Since the involution exchanges 
the two cases, we can assume that u corresponds to the first case. Knowing the term A u , it is 
possible to compute the contributions of oj and s(u) to all of the edges adjacent to v. Indeed, 

• The probability of s(oj) is equal to x(p)y2 times the probability of u (due to the fact that 
there is one additional open edge and one additional loop); 

• Windings of the curve can be expressed using the winding at A. For instance, the winding at 
B in the configuration u is equal to the winding at A minus a tt/2 turn. 

The contributions are given as: 



configuration 


A 


C 


B 


D 


OJ 


Auj 


e i7r / 2 ^ 


e- i7r / 4 A, 




s(co) 


x(p)V2A 0J 








e i7r / 4 x(p)\/2A, 



Using the identity e in/A + e _i7r/4 
edges around v. 



\/2, we deduce (2.5) by summing the contributions of all the 

□ 



The previous lemma provides us with one linear relation between values of F for every vertex 
inside the domain. However, there are approximately twice as many edges than vertices in G 
so that these relations do not completely determine the value of F. The next lemma is therefore 
crucial since it decreases the number of possible values for F; roughly speaking, it states that the 
complex argument (modulo n) of F(e) is determined by the orientation of the edge. 



Lemma 2.4. F(e) belongs to M. (resp. e~ m ^M., MR or e m ^M.) on edges e pointing in the same 
direction as the ending edge e& (resp. edges pointing in a direction which forms an angle tt/2, it and 
3n/2 with e^). 

Proof. The winding at an (oriented) edge can only take its value in the set Wq + 2ttZ where Wo 
is the winding at e of an arbitrary possible interface passing through e. Therefore, the winding 
weight involved in the definition of F is always proportional to e lW °/ 2 with a real-valued coefficient, 
and thus the complex argument of F is equal to Wo/2 or Wo/2 + n. Since Wo is exactly the angle 
between the direction of e and that of eft, we obtain the result. □ 
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The observable in strips. The definition of F can be extended to the case of the strip. Indeed, 
the loop representation extends in this setting; the c/>^' _00 -probability of having an infinite cluster 
is 0: for fixed £, the model is essentially one-dimensional and it is a simple exercise to prove that 
it must be subcritical. Hence, there is a unique interface going from +00 to —00, which we call 7. 
We define 

F(e) := ^ 

where W 7 (e, —00) is the winding of the curve between e and —00. This winding is well-defined up 
to an additive constant, and we set it to be equal to for edges of the bottom side which point 
inside the domain. It is easy to see that F is the limit of observables in finite boxes, so that the 
properties of fermionic observables in Dobrushin domains carry over to the infinite- volume case. In 
particular, the conclusions of the previous three lemmas apply to it as well. 
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3 Proof of Theorem [T] 

The proof consists of three steps: 



We first prove using Lemmas 2.3 and 2.4 that the observable decays exponentially fast when 



P < Psd in a we ll chosen Dobrushin domain (namely a strip with free boundary condition 



on the top and wired boundary condition on the bottom). Lemma 2.2 then implies that the 
probability that a point on the top of the strip is connected to the bottom decays exponentially 
fast in the height of the strip. 

• We derive exponential decay of the connectivity function for the infinite- volume measure with 
free boundary conditions from the first part. 

• Finally, we show that exponential decay implies that the random-cluster model is subcritical 
when p < p s d, and that its dual is supercritical. This last step concludes the proof of 
Theorem Q] and is classical. 

In the proof, points are identified with their complex coordinates. 



Step 1: Exponential decay in the strip. Let p < p s d and consider the random-cluster model 
on the strip Si of height £ > with wired boundary condition on the bottom and free boundary 
condition on the top. Define e k and e^+i to be the north-west-pointing sides of the diamonds 
associated to the points ik and i(k + 1), respectively. Label some of the edges around these two 
diamonds as x, x', x", y and y 1 as shown in Figure |4| 



Lemmas 2.3 and 2.4 have a very important consequence: around a vertex v, the value of the 



observable on one edge can be expressed in terms of its values on only two other edges. This can 



be done by seeing the relation given by Lemma 2.3 as a linear relation between four vectors in the 
plane M 2 , and applying an orthogonal projection to a line orthogonal to one of them (which can be 



chosen using Lemma 2.4). One then gets a linear relation between three real numbers, but using 
Lemma |2.4| "in reverse" shows that this is enough to determine any of the corresponding three 
(complex) values of the observable given the other two. 

For instance, we can project ( |2.4[ ) around V\ orthogonally to F(y), so that we obtain a rela- 
tion between projections of F(x), F{x') and F(ei s +i)- Moreover, we know the complex argument 
(modulo 7r) of F for each edge so that the relation between projections can be written as a relation 
between F(x), F(x') and F^e^+i) themselves. This leads to 

e-' m/4 F(x) = cos(tt/4 - a)F(e k+1 ) - cos(tt/4 + a)e" i7r/2 'F{x'). (3.1) 
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Figure 4: Left: The labelling of edges around used in Step 1. Right: A dual circuit 
surrounding an open path in the box [— 02, a2] 2 . Conditioning on to the most exterior such 
circuit gives no information on the state of the edges inside it. 



Applying the same reasoning around t>2, we obtain 

e~ in/i F(x) = cos(^/4 + a)F(e k ) - cos(tt/4 - a)e- in/2 F(x"). 
The translation invariance of implies 

F(x') = F(x"). 

Moreover, symmetry with respect to the imaginary axis implies that 

F(x) = e^^Fjx 7 ) = e- i7r/i F(x'). 



(3.2) 



(3.3) 



(3.4) 



Indeed, if for a configuration to, x belongs to 7 and the winding is equal to W, in the reflected 
configuration to', x' belongs to 7(0/) and the winding is equal to tt/2 — W. 
Plugging (3.3) and (3.4) into ( |3.1[ ) and ( |3.2[ ), we obtain 

_ i7r/4 1 + cos(tt/4 + a) [1 + cos(tt/4 + a)} cos(tt/4 + a) 

K k+l) cos(vr/4-a) { ' [1 + cos(tt/4 - a)] cos(vr/4 - a) [ kh 

Remember that a{p) > since p < p s d, so that the multiplicative constant is less than 1. Using 
Lemma 2.2 and the previous equality inductively, we find that there exists c\ = c\{p) < 1 such 
that, for every t > 0, 

A = \F(e e )\ = c{\F( ei )\ < c{ , 
where the last inequality is due to the fact that the observable has complex modulus less than 1. 



Step 2: Exponential decay for (p^ when p < p^. Fix again p < p S( i- Let and recall 

that (fi® N := (f)® 2 r/vrjv] 2 converges to the infinite- volume measure with free boundary conditions 
(ftp when N goes to infinity. 

Consider a configuration in the box [— N, N] 2 , and let A mSLX be the site of the cluster of the 
origin which maximizes the £°°-norm max{|xi|, |x2|} (it could be equal to N). If there is more 
than one such site, we consider the greatest on6 in lexicographical order. Assume that -Amax 

equals 

a = a\ + \d2 with a2 ^ I ai| (the other cases can be treated the same way by symmetry, using the 
rotational invariance of the lattice) . 

By definition, if A miiX equals a, a is connected to in [—02, o^] 2 . In addition to this, because of 
our choice of the free boundary condition, there exists a dual circuit starting from o + i/2 in the 
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dual of [— d2, 02] 2 (which is the same as L* n [—02 — 1/2, 02 + 1/2] 2 ) and surrounding both a and 0. 
Let r be the outermost such dual circuit: we get 



$U0W = «) = 5>lU( a °l r = ^X^v( r = T)> (3-5) 

7 

where the sum is over contours 7 in the dual of [— a2;«2] 2 that surround both a and 0. 

The event {r = 7} is measurable in terms of edges outside or on 7. In addition, conditioning 
on this event implies that the edges of 7 are dual-open. Therefore, from the domain Markov 
property, the conditional distribution of the configuration inside 7 is a random-cluster model with 
free boundary condition. Comparison between boundary conditions implies that the probability of 
{a 0} conditionally on {r = 7} is smaller than the probability of {a f-> 0} in the strip S a2 with 
free boundary condition on the top and wired boundary condition on the bottom. Hence, for any 
such 7, we get 

4>l N {a o o|r = 7) < C 2 "°> *+ 0) = C;°> z ) < C T = 4 a|/2 

(observe that for the second measure, Z is wired, so that {a 0} and {a o Z} have the same 



probability). Plugging this into (3.5), we obtain 



<^ max = a) < ^(r = 7) < 4 a|/2 . 

7 

Fix n ^ N. Since ci < 1, we deduce from the previous inequality that there exist two constants 
< C2, C2 < 00 such that 

^ j7V (0oZ 2 \[-n,n] 2 ) < J2 <fiA A *»* = a )< E c i°' /2 < C2e- C2n 

ae[-7V,iV] 2 \[-n,n] 2 a^[-n,n] 2 

Since the estimate is uniform in N, we deduce that 

<^(0 o Z 2 \ [-n, n] 2 ) ^ C 2 e- C2n . (3.6) 



Step 3: Exploiting exponential decay. The inequality p c ^ p S( i follows from (3.6) since 
exponential decay prevents the existence of an infinite cluster for cj)® when p < p S( i- 

In order to prove that p c ^ p s d, we use the following standard reasoning. Let A n be the event 
that the point (n, 0) is in an open circuit which surrounds the origin. Notice that this event is 



included in the event that the point (n, 0) is in a cluster of radius larger than n. For p < p s d, (3.6) 
implies that the probability of A n decays exponentially fast. The Borel-Cantelli lemma shows that 
there is almost surely only a finite number of values of n such that A n occurs. In other words, there 
is only a finite number of open circuits surrounding the origin, which enforces the existence of an 
infinite dual cluster. It means that the dual model is supercritical whenever p < p s d- Equivalently, 
the primal model is supercritical whenever p > p sc i, which implies p c ^ p s d- □ 



4 Proof of Theorem [2] 

In this section, we compute the correlation length in all directions. In [14J, Messikh noticed that 
this correlation length was connected to large deviations for random walks and asked whether there 
exists a direct proof of the correspondence. Indeed, large deviations results are easy to obtain for 
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random walks, so that one could deduce Theorem [2] easily. In the following, we exhibit what we 
believe to be the first direct proof of this result. 

An equivalent way to deal with large deviations of the simple random walk is to study the 
massive Green function G m , defined in the bulk as 



G m (x,y) :=E*[j2™ n tx n . 



where K x is the law of a simple random walk starting at x. 

The correlation length of the two-dimensional Ising model is the same as the correlation length 
for its random-cluster representation so that we will state the result in terms of the random-cluster. 
We use the parameters p and a = a(p) without revealing the connection with f3 in the notation. 

Proposition 4.1. For p < p s d and any a G L, 

- lim - log 0°(O O na) = - lim - log G m (0, raa) (4.1) 

n— >oo n ' n— >oo n 



where m = cos[2a(p)] - the value of a(p) is given by (2.3). 



In p3], the statement involves Laplace transforms but we can translate it into the previous 
terms. Moreover, the mass is expressed in terms of /3, but it is elementary to compute it in terms 
of q. Theorem [2] follows from this proposition by first relating the two-point functions of the Ising 
and q = 2 random-cluster models as was mentioned earlier, and then deriving the asymptotics 
of the massive Green function explicitly — the details can be found for instance in the proof of 
Proposition 8 in |14j . 

Before delving into the actual proof, here is a short outline of the strategy we employ. We have 
already seen exponential decay in the strip, which was an essentially one- dimensional computation; 
we want to refine it into a two-dimensional version for correlations between two points and a in 
the bulk, and once again we use the observable to estimate them. The basic step, namely obtaining 
local linear relations between the values of the observable, is the same, although it is complicated 
by the lack of translation invariance. The point is that the observable is massive harmonic when 



P 7^ Psd (see Lemma 4.2 below). Since G m (-,-) is massive harmonic in both variables away from 
the diagonal x = y, it is possible to compare both quantities. 

The main problem is that we are interested in correlations in the bulk. The observable can be 
defined directly in the bulk (see below) but it provides only a lower bound on the correlations. In 
order to obtain an upper bound, we have to introduce an "artificial" domain (that will be T(a) 
below), which needs two features: the observable in it can be well estimated, and at the same 
time correlations inside it have comparable probabilities to correlations in the bulk. For the second 
one, it is equivalent to impose that the Wulff shape centered at and having a on its boundary is 
contained in the domain in the neighborhood of a; from convexity, it is then natural to construct 
T(a) as the whole plane minus two wedges, one with vertex at and the other with vertex at a. 

The proof is rather technical since we need to deal with the behavior of the observable on the 
boundary of the domains. This was also an issue in Smirnov's proof. At criticality, the difficulty 
was overcome by working with the discrete primitive H of F 2 . Unfortunately, there is no nice 
equivalent of H to work with away from criticality. The solution is to use a representation of F 
in terms of a massive random walk. This representation extends to the boundary and allows to 
control the behavior of F everywhere. 
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Proof. Let p < p S( [. Without loss of generality, we can consider a = (01,02) £ L satisfying 
02 ^ 01 ^ 0. In the proof, we identify a site u of L with the unique side e u of the associated black 
diamond which points north-west. In other words F(u) and {u G 7} should be understood as F(e u ) 
and {e n S 7} — notice that this differs from the notation used in [16]. 



The lower bound. Consider the observable F in the bulk defined as follows: for every edge e 
not equal to eo, 

F(e):=We^ e < e °)l ee7 V (4.2) 



where 7 is the unique loop passing through eo- Note that this definition is justified by the fact that 
p is subcritical, and that it immediately implies that 



$(0*+a) > \F(a)\. 



(4.3) 



We mention that F is not well defined at eo- Indeed, eo can be thought of as the start of the loop 
7 or its end. In other words, F is multi- valued at eo, with value 1 or -1. 



Lemma 2.3 can be extended to this context following a very similar proof, but taking into 
account that F is multi-valued at eo- More precisely, let eo = xy. Around any vertex v ^ {x,y} 



the relation in Lemma 2.3 still holds; besides, 



(F(SE) + 1 = e iQ ^ [F(SW) + F(NE)} if v = y 
\F(SW) + F(NE) = e iQ ( p ) [-1 + F(SE)] if v = x 

where the NE (resp. SE, SW) is the edge at v pointing to the north-east (resp. south-east, south- 



west). In other words, the statement of Lemma 2.3 still formally holds if we choose the convention 
that F{e$) = 1 when considering the relation around x, and F{eo) = —1 when considering the 
relation around y. 



One can see that Lemma 2.4 is still valid. In fact, the two lemmas imply that F is massive 
harmonic: 

Lemma 4.2. Let p < p s d and consider the observable F in the bulk. For any site X not equal to 
0, we have 

A a F(X) := -^[ F (W)+F(S) + F(E) + F(N)]-F(X) = 0, 
where W , S, E and N are the four neighbors of X. 

Proof. Consider a site X inside the domain and recall that we identify X with the corresponding 
edge of the medial lattice pointing north-west. Index the edges around X in the same way as in 
Case 1 of Figure [5| By considering the six equations corresponding to vertices that end one of 
the edg (being careful to identify the edges A, B, C and D correctly for each of the 

vertices), we obtain the following linear system: 



F{X) + F(yi) 
F{y 2 ) + F{x x ) 
F(S) + F(x 2 ) 
F{x 3 ) + F(x 4 ) 
F(E) + F(x 5 ) 
F(x 6 ) + F{y 6 ) 



[F( Xl 
[F(x 2 

[F(V3 
[F(y 4 
[F(x 4 
[F(x 5 



+ F(x 6 )] 
+ F(W)] 
+ F(x 3 )] 
+ F(X)] 
+ F(y 5 )} 
+ F(N)] 
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case 1 



case 2 



case 3 



case 4 



Figure 5: Indexation of the edges around vertices in the different cases. 



Recall that by definition, F(X) is real. For an edge e, denote by /(e) the projection of F(e) on 
the line directed by its argument (R, e 17r//4 M, iM. and e -17r//4 IR). By projecting orthogonally to the 
F(yi), i = 1 . . . 6, the system becomes: 



f fix) 


= cos(7r/4 + 


a)f(x 1 ) 


+ 


cos(7r/4 — a)f(xe) 


(1) 




= cos(7r/4 + 


a)f(x 2 ) 


+ 


cos(^/4 - a)f (W) 


(2) 


/&») 


= cos(7r/4 — 


<*)f(S) 


+ 


- cos(7r/4 + a)f(x 2 ) 


(3) 


/(*) 


= cos(7r/4 + 


a)f(x 3 ) 


+ 


cos(7r/4 — a)f(x4) 


(4) 




= cos(7r/4 + 


a)f(E) 


+ 


cos(7r/4 — a)f{x§) 


(5) 


I /(a*) 


= — cos(7r/4 


- at)f(x 5 ) 


+ 


cos(^/4 + a)f{N) 


(6) 



By adding (2) to (3), (5) to (6) and (1) to (4), we find 



f(x 3 ) + f(x 1 ) = 
f(x 6 ) + f(x 4 ) = 
2/PO = 



cos(n/A- a)[f(W) + f(S)\ 
cos(7r/4 + a)[/(^) + /(JV)] 

cos(vr/4 + a)[f(x 3 ) + f{x{)} + cos(vr/4 - a)[/(x 6 ) + /(x 4 )] 



(7) 
(8) 
(9) 



Plugging (7) and (8) into (9), we obtain 

2f(X) = cos(^/4 + a) cos(tt/4 - a)[f(W) + f(S) + f(E) + f(N)}. 

The edges X, . . . , N are pointing in the same direction so the previous equality becomes an equality 
with F in place of / (use Lemma 2.4). A simple trigonometric identity then leads to the claim. □ 



Define the Markov process with generator A a , which one can see either as a branching process 
or as the random walk of a massive particle. We choose the latter interpretation and write this 
process (X n , m n ) where X n is a random walk with jump probabilities defined in terms of A a — the 
proportionality between jump probabilities is the same as the proportionality between coefficients 
- and m n is the mass associated to this random walk. The law of the random walk starting at x 
is denoted ¥ x . Note that the mass of the walk decays by a factor cos 2a at each step. 

Denote by r the hitting time of 0. The last lemma translates into the following formula for any 
a and any t, 

F(a)=E a [F(X tAr )m tAT }. (4.4) 



The sequence (F(Xt)mt)t^r is obviously uniformly integrable, so that (4.4) can be improved to 

F{a) =E a [F(X T )m T }. (4.5) 
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Equations (4.3), (4.5) together with Lemma 4.3 below give 



0°(O <-> a) ^ T^Gfcoete^a), 
|a| 

which implies the lower bound. 

Lemma 4.3. There exists c > such that, for every a in the upper-right quadrant, 

\E a [F(X T )m T }\ > ^G cos2Q (0,a). 
|a| 

Proof. Recall that F(X T ) is equal to 1 or -1 depending on the last step the walk takes before 
reaching 0. Let us rewrite E a [F(X T )m T ] as 

E a [m T 1 {Xt _ 1= W or S}] - 1 {Xt _ 1= N or E}]- 

Now, let A a be the line y = —x, and let T be the time of the last visit of A Q by the walk before 
time t (set T = oo if it does not exist). On the event that X T _i = W or S, this time is finite, 
and reflecting the part of the path between T and r across A a produces a path from a to with 
X T -i = E or N. This transformation is one-to-one, so summing over all paths, we obtain 



HX T . 1= W or S}} ~ 1{X T _ 1= ^V or E}] 



E a \m T 1 



{X r _i=7V or 



£}l{T=oo}] 



which in turn is equal to —E a [m T l^ =00 y]. General arguments of large deviation theory imply that 
E a [m T l|j' =00 }] > AG COS 2a(0, a) for some universal constant c. □ 



case 3: vertical part of L + (w) 



case 4: site w 




AAA 



case 2: horizontal part of L + (w) 
case 1: interior of T(w) 



wired arc L of T(w) 



Figure 6: The set T(w). The different cases listed in the definition of the Laplacian are 
pictured. 



The upper bound. Assume that is connected to a in the bulk. We first show how to reduce 
the problem to estimations of correlations for points on the boundary of a domain. 

For every u = u\ + fu2 and v = v\ + if 2 two sites of L, write u -< v if u\ < v\ and U2 < V2- This 
relation is a partial ordering of L. We consider the following sets 

L+(u) = {x G L : u -< x} and L~ = {x G L : x -< 0}; 
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and 

T{u) = L\ (L+(u) UL"). 

In the following, L + {u) and L~ will denote the interior boundaries of T{u) near L + («) and L~ 
respectively, see Figure [6| The measure with wired boundary conditions on L~ and free boundary 
conditions on h + (u) is denoted </>t(u)- 

Assume that a is connected to in the bulk. By conditioning on w which maximizes the partial 
^-ordering in the cluster of (it is the same reasoning as in Section 3), we obtain the following: 

(jpja <H> 0) ^ \^ 4>t(w)(w ^ Cz\a\ max i/i TM (wf)>r) (4.6) 

' «)!-a,|ui|sgC3|a| v ; 

for C3, C3 large enough. The existence of C3 is given by the fact that the two-point function decays 
exponentially fast: a priori estimates on the correlation length show that the maximum above 
cannot be reached at any w which is much further away from the origin than o, and even that 
the sum of the corresponding probabilities is actually of a smaller order than the remaining terms. 



Summarizing, it is sufficient to estimate the probability of the right-hand side of (4.6). 



Observe that w is on the free arc of T(w), so that, harnessing Lemma 2.2 we find 

<t> T{w) {w ^ L-) = \F(w)\, (4.7) 
where F is the observable in the Dobrushin domain T(w) (the winding is fixed in such a way that 



it equals at e w ). Now, similarly to Lemma 4.2, F satisfies local relations in the domain T(w): 



Lemma 4.4. The observable F satisfies A a F = for every site not on the wired arc, where 
the on T(w) is defined by the following relations: for all g : T(w) h-> E, 

(g + A a g)(X) is equal to: 

cos 2o^ 

— ^ — [5(^0 + 9(8) + + g(N)] inside the domain; 

cos2a r . . cos(7r/4 + a) , , , . r j., \ 

[g{W) + g(S)] + v ' — —;g(E) on the horizontal part ofL + (w); 



2(1 + cos(tt/4 - a)) L v ' y n l + cos(7r/4-a 

cos2a , .„,. cos(7r/4 + a) ,_ T . , . , rr j., \ 

wfr- r^- A yr\g (W )+g{S)} + , . g(N) on the vertical part of L + {w); 

2(1 + cos(7T/4 — a)) 1 + cos(7r/4 — a) 

cos 2a r cos(7r/4 — a) r /rnX , , TX1 
-^ r [g(W)+g(S)) + ^^ 2 >-[g{E) + g{N)] atw, 

with N, E, S and W being the four neighbors of X. 

Proof. When the site is inside the domain, the proof is the same as in Lemma |4.2| For boundary 
sites, a similar computation can be done. For instance, consider Case 2 in Fig. [5j Equations (3) 



and (7) in the proof of Lemma 4.2 are preserved. Furthermore, Lemma 2.2 implies that 



f(X) = f(x l ) = <j> T{w) (X^L-) 



and similarly f{x^) = f(E) (where / is still as defined in the proof of Lemma 4.2). Plugging all 



these equations together, we obtain the second equality. The other cases are handled similarly. □ 

Now, we aim to use a representation with massive random walks similar to the proof of the 
lower bound. One technical point is the fact that the mass at w is larger than 1. This could a 
priori prevent {F{Xt)mt)t from being uniformly integrable. Therefore, we need to deal with the 
behavior at w separately. Denote by n the hitting time (for t > 0) of w, and by r the hitting time 
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of L . Since the masses are smaller than 1, excepted at w, {F{Xt)mt)t< T hT 1 is uniformly integrable 
and we can applying the stopping theorem to obtain: 

F(w) = E w [F(X TAT1 )m rAT1 ] = E w [F(X T1 )m T1 t T1<T ]+E w [F(X T )m T t T<T1 ]. 

Since X Tl = w, the previous formula can be rewritten as 

E w [F(X T )m T l T<T1 ] 



F{w) 



l-E«>(m ri l T1<T ) ' 



(4.8) 



When w goes to infinity in a prescribed direction, [1 — E w (m Tl t T1<T )] converges to the analytic 
function h : [0, 1] — > M, p H > 1 — E w (m T1 ) (since the function is translation-invariant). The function 
h is not equal to when p = 0, implying that it is equal to for a discrete set V of points. In 
particular, for p ^ V, the first term in the right hand side stays bounded when w goes to infinity. 
Denoted by C4 = C±(jp) such a bound. Recalling that \F\ ^ 1 and that the mass is smaller than 1 
except at w, (4.8) becomes 



\F(w)\ < C 4 \E w [F(X T )m T t T<T1 ]\ < E w [m T l T<Tl ] 

< C 4 ^ E x [(cos2a) T t T<T1 l {(Xt) avoids i+H} ] ^C 4 J2 G 



cos 2a 



(0,x) 



(4.9) 
(4.10) 



w~<x w^,x 

where the last inequality is due to the fact that we release the condition on avoiding L + (w). 
Finally, it only remains to bound the right hand side. From (4.10), we deduce 

\F{w)\ < C 5 \w\G cos2a (0,w) (4.11) 

where the existence of C5 is due to the exponential decay of G COS 2o(-, •) and the fact that G COS 2a(0, x) ^ 
Cr C os2a(0, w ) whenever w <x. We deduce from (4.6), (4.7) and (4.11) that 



P (O o a) < C 3 C 5 \a\ 



max 



■1TC 



G m (0,w) < C 6 |a| 2 G m (0,a) 



(4.12) 



Taking the logarithm, we obtain the claim for all p < p S( i not in the discrete set V. The result 
follows for every p using the fact that the correlation length is increasing in p. □ 
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